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Abstract 



We introduce a momentum subtraction scheme which obeys the power 
counting of Kaplan, Savage, and Wise (KSW), developed for systems with 
large scattering lengths, a. Unlike the power divergence subtraction scheme, 
coupling constants in this scheme obey the KSW scaling for all /j,r > 1/a. 
We comment on the low-energy theorems derived by Cohen and Hansen. We 
conclude that there is no obstruction to using perturbative pions for momenta 
p > m n . 
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Effective field theory is a useful method for describing two-nucleon systems |J . Recently, 
Kaplan, Savage, and Wise (KSW) devised a power counting that accounts for the effect of 
large scattering lengths. With this power counting the dimension six four-nucleon operators 
are treated non-perturbatively, while pion exchange as well as higher dimension operators 
are treated perturbatively. A key feature of the KSW approach is the use of the power 
divergence subtraction (PDS) renormalization scheme. PDS gives the coefficients of contact 
interactions power law dependence on the renormalization scale, which makes the KSW 
power counting manifest. For > 300 MeV, the power counting in PDS is no longer 
manifest @, and KSW conclude that application of the theory is restricted to momenta less 
than 300 MeV. 

In Ref. 0, it is emphasized that the power counting for large scattering lengths should be 
scheme independent. Here we introduce a momentum subtraction scheme that is compatible 
with the KSW power counting. A similar scheme is applied to the theory without pions in 
Ref. |4]]. This scheme will be referred to as the OS scheme, since in a relativistic field theory 
it would be called an off-shell momentum subtraction scheme. One attractive feature of 
the OS scheme is that the KSW power counting is manifest for all fiR > 1/a, where a is 
the scattering length. Since the breakdown of the power counting is an artifact of the PDS 
scheme, the range of the effective field theory with perturbative pions may not be limited 
to 300 MeV. 

Another attractive feature of the OS scheme is that the amplitudes are manifestly fiR 
independent at each order in the expansion. In PDS, any fixed order calculation has residual 
fin dependence which is cancelled by higher order terms. The rather strong [ir dependence 
exhibited by next-to-leading order calculations of the phase shifts led authors [[||| to con- 
clude that PDS results have to be fine tuned to fit the data. In fact, it is possible to fix this 
problem with PDS by modifying the treatment of Co(/zr). 

We begin our discussion by recalling the Lagrangian with pions and nucleons @ 

£ n = £ Tr d^) + ^ Tr(m 9 £ + m q tf) - ^ N^a^d^ - 
o 4 Z 

+ rf(iD + ^)N- C^(N t P^N)\N t P^N) (1) 

- D { 2 s) ujTi(m^(N T P} s) Ny(N T P} s) N) + .... 

Here qa = 1.25 is the nucleon axial- vector coupling, / = 130 MeV is the pion decay constant, 
= \{im q ^ + ^m q ^), where m q = diag(m M ,mrf) is the quark mass matrix, and = 
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(N T P} s) Ny(N T Pl s) V 2 iV) + h.c. 



w(m u + rrid)- The matrices P^ s project onto states of definite spin and isospin, and the 
superscript s denotes the partial wave amplitude mediated by the operator. This paper will 
be concerned only with S-wave scattering, so s = 1 Sq, 3 Si. The ellipsis in Eq. (P denotes 
terms that are higher order in the expansion. 

The four nucleon couplings Co, C 2 , and D 2 in Eq. (JJ) are bare parameters. In general, 
a bare coupling, C bare , can be separated into a renormalized coupling, C(//r), and counter- 
terms: 

oo 

c bare = c finite - 5 uv c , c flnite = c(n R ) - J2 $ n c(m) • (2) 

n=0 

The counterterms have been divided into two classes. The first, which have the superscript 
uv, contain all ultraviolet divergences and are \Lr independent. Defining a renormalization 
scheme amounts to making a choice for the finite counterterms, S n C(fiR,). The superscript 
n indicates that 5 n C is included at tree level for a graph with n loops. For example, at 
two loops, we have two loop diagrams with renormalized couplings at the vertices, one loop 
diagrams with a single S l C counterterm, and a tree level diagram with 5 2 C. Let Q denote 
a typical momentum characterizing the process under consideration. For nucleon-nucleon 
scattering we take p ~ Q and m n ~ Q, where p is the center-of-mass momentum. The 
theory is an expansion in Q/A where A is the range of the effective field theory. Taking 
fJ'R ~ Qi vertices with C (//r) scale as 1/Q, while vertices with C 2 p 2 or D 2 m 2 scale as Q°. A 
typical loop gives one power of Q, so Cq(hr) vertices are included to all orders. This sums 
all corrections that scale as (Qa) n [[|. Note that since the pion has been included explicitly 
in the Lagrangian we expect that the scale of short distance physics, A, should not be set 
by m n , but by higher mass resonances which have not been included in the theory. 

PDS is one scheme in which the KSW power counting is manifest. In PDS, we first let 
d = 4 — 2e and define the counterterms 5 UV C to subtract 1/e poles. As in the MS scheme, the 
dimensional regularization parameter /i is set to //r. Next one takes d = 3 and defines the 
finite counterterms, 5 n C(/j,R) to subtract the l/(d— 3) poles in the amplitude. Graphs which 
renormalize a given coupling are those whose vertices have the right number of derivatives 
and powers of m 2 . When calculating 5 n C (fifi) and 5 n C 2 (/iR), we can take m n = since 
counterterms proportional to m 2 renormalize coefficients like D 2 {^r)- After making these 
subtractions, the amplitude is continued back to four dimensions, so d — 3 — > 1. 

In the 1 Sq channel, exact expressions for the PDS beta functions can be obtained [El. 
For C ([ir) we have 
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where A^ty = ^f 2 /{Mg\) = 300 MeV. if is a constant fixed by the boundary condition, 
and choosing C (0) = Ana/M gives if = l/(a + 1/A^n)- We see that the scaling for Co(/xr) 
changes for ~ 300 MeV 0. A simple shift, C (/i#) — > C (/xr) — g\/2f 2 , results in a 
coupling that scales as 1/fiR for all fiR> 1/a. Physically, this shift corresponds to summing 
the short distance (m^ = 0) contributions from potential pion exchange to all orders. For 
the 3 5i channel, this summation is not possible because of ultraviolet divergences of the form 
p 2 /e and, in fact, there are unknown corrections to the 3 Si beta functions at each order in 
Q. This is demonstrated in Ref. J7| by explicit computation of the counterterms. When the 
beta function is not exactly known, the large fiR behavior is ambiguous. For example, the 
PDS beta function for Co(jUr) is 



dC (fi R ) Mfi R 



,2 



dfi R An 

Two solutions which satisfy this equation to order Q° are 

4:71 



+ 0(Q) . (4) 



Coins) 



2a\iR - 2fi g - J I + Afi 2 - 4// ff (l - afi R ) 



M 2fi R {\ - afiR - fig) 



CM = -% - r — -J— ■ — — , (5) 



Pi- 



1 + 2/(aA NN ) exp (— 2/i 



where fi g = fiR/A NN and we have chosen C (0) = Ana/M. The first solution is obtained 
by computing the counterterms 5 n Co(fiR) to order Q° and summing them. This solution 
falls as 1/fiR for all fiR> 1/a, and is numerically close to the — ► solution. The second 
solution is obtained by truncating and solving Eq. M). This solution approaches a constant 
as fiR —>■ oo. The two solutions both solve the beta function to order Q° but have very 
different large fi R behavior. In the OS scheme, there is no ambiguity since at a given order 
in Q the running of all the coupling constants that enter is known exactly. 

In the OS momentum subtraction scheme, the renormalized couplings are defined by 
relating them to the amplitude evaluated at the unphysical momentum p = ifiR. We start 
by dividing up the full amplitude as A = J2m=o A^" 1 " 1 ^. Here A^ m ~^ contains the Feynman 
diagrams that scale as Q™ 1 " 1 . As in PDS, we can take — > when computing the 
counterterms for Co(fiR) and C^/xr). S uv C2 m is first defined to subtract all four dimensional 
1/e poles. The definition for the renormalized coupling is then 
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FIG. 1. Renormalization conditions for Co(/xr) and C 2 ([ir) in the OS scheme, ivl* 1 ) is the 
four point function with l7o(/xr) and 5 n Co(nn) vertices, evaluated between incoming and outgoing 
1 Sq or 3 Si states. The amplitude contains graphs with one C 2 or one potential pion dressed 
with Cq bubbles. 
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(6) 



This condition is to be imposed order by order in the loop expansion so that the graphs at 
n loops determine 5 n C (fiji). For instance, the couplings C (/iR) and C 2 (/i#) are defined by 
the renormalization condition in Fig. [I]. Summing the counterterms, we find 
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Although it may seem that the piece of C 2 (/jr) that goes as will spoil the power 

counting for low momentum, in fact, the part dominates entirely for /xr > 1/a, since 



finite _ /^(finite „2 
~ a, G 2 ~ d • 



In the OS scheme, D 2 {^lr) will be calculated as follows. First m^/e poles are subtracted. 
The renormalized coupling is then defined by 



iA(D 2 



-i D 2 (fJ, R )ml 



(8) 



where A(D 2 ) contains terms in the amplitude that are proportional to m%, as well as analytic 
in m%. Thus, at Q° graphs with a single D 2 {hr) or potential pion and any number of Co 
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vertices contribute to A(D 2 ). Note that we only keep terms that are analytic in m 2 because 
it seems unnatural to put long-distance nonanalytic contributions that come from pion 
exchange into the definition of the short distance coupling D 2 (hr) 0. For example, one 
potential pion exchange gives a (m 2 /p 2 )ln(l + Ap 2 /m 2 ) term. Putting this into D 2 {^r) 
would give it both a branch cut at fiR = as well as explicit dependence on the scale 

m n , therefore this term is not included in A(D 2 ). Solving for the counterterms and summing 
we find 

D 2 (jj,r) 
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Co^r) 2 (C^) 2 sAsttP " 
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where /z 2 = /i 2 , exp ( 1 



647r 2 / 2 Df nitc 



(9) 



M 2 g 2 A (C$ nite ) 2 

With ~ Q ~ fiR, D 2 (fiu)ml ~ Q°. The scale ft must be determined by fitting to data. 
The PDS solution for D 2 ([ir) does not have the —1 in square brackets. 

We now compare the amplitudes in the PDS and OS schemes. At order Q°, the ampli- 
tudes in the ^So and 3 S± channels have the same functional form. The amplitudes in PDS 
are calculated to this order in Ref. H. We find 
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(10) 
(11) 



(12) 
(13) 



Note that the last term in Eq. (^) has a factor of 1 /2 instead of a 1 as in M since we have 
made a different finite subtraction. The terms in braces are long distance pion contributions, 
and are the same in the PDS and OS schemes. By substituting Eqs. (|?|) and (Jfty into the 
OS amplitude one can verify that it is fiR independent. In contrast, in PDS the amplitude 
is only [ir independent to order Q°. 
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To obtain a good fit to the scattering data at low momenta, two constraints must be 
approximately satisfied: 

1 _ a® 

5 r 

a 



lim pcot5(p) 
p— >o 



[A(-i)] 1 



0, 



(14) 



where 7 = MC f^j The second constraint ensures that the next-to-leading order ampli- 

tude does not have a spurious double pole. In the PDS scheme, is not /i# independent, 
so the extracted parameters can not simultaneously satisfy the renormalization group equa- 
tions and give a good fit. This is the origin of the large dependence on observed in 
Ref. ||. In PDS the second constraint gives 

-M\ 2 rl 1, f u%\ 2 / 4?r 
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4tt 



(15) 



MC (fij 



while in OS we find 

-D 2 (ij, r ) 







Co (A* 



2f 2 UtJ 



2 



47T 



(16) 



2"Vm 2 7 m^MCoW 
where terms of order 7 2 /m 2 have been dropped. With the convention for D 2 ([ir) in Ref. 
the first term in square brackets in Eq. ([15]) is a 1, so they find that D^m^) is close to zero. 
The last term in Eq. (|15|) induces the large hr dependence. On the other hand Eq. ( |I~6"D is 
\i r independent. 

Fits are performed to the Nijmegen phase shift || data between 7 and 100 MeV for both 
the 1 5'o and 3 S\ channels. The phase shifts have an expansion of the form 5 = 5^ + 5^ + 
e>(Q 2 /A 2 ), where § 
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4tt 1 + iftA(-V 



(17) 



The fits were weighted towards low momentum since the theoretical error is smallest there. 
The results are shown in Fig. |2|. We chose to fit using the amplitudes in the OS scheme, but 
an equally good fit is obtained in PDS. The parameters C ([1r), C 2 ([j,r), and D 2 (fi R ) were 
extracted for different values of fiR in both schemes, and the conditions in Eq. (p~5| ) and ( |16| ) 
were found to be well satisfied. For instance, taking fiR = m n = 137 MeV in the OS scheme 
we find 
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FIG. 2. Fit to the phase shift data weighted toward low momentum. The solid line is the 
Nijmegen fit to the data ||, the long dashed line is the order 1/Q result, and the short dashed line 
is the order Q° result. 

The fit gives scattering lengths a( 1 S'o) = — 23.3fm, and a( 3 Si) = 5.41 fm. Plugging the 
parameters extracted from the fit into the right hand side of Eq. fll6|) gives ~ 10~ 2 for 
both channels. In the OS scheme, the renormalization group equations are obeyed to a few 
percent accuracy since the amplitude is explicitly hr independent. 

It is important to realize that the fits do not unambiguously determine the values of 
Cq(hr) and D 2 (fiR). The coefficient of the four nucleon operator with no derivatives is 
C 6ar£ ' J r m l x D\ are . Once we switch to renormalized coefficients, it is not clear how to divide the 
couplings into a nonperturbative and perturbative piece. In Ref. 0, Co(aO is summed to all 
orders, while other authors |5]|| treat both coefficients non-perturbatively. In fact, in order 
to do a chiral expansion, m 2 r D 2 (^_R) should be treated perturbatively. Since m 2 r D 2 (/^i?) ~ Q°, 
this is consistent with the power counting and the renormalization group equation. 

On the other hand, there is some freedom in dividing Cq(hr) into nonperturbative and 
perturbative pieces: Cq(iir) = Cq P {^r) + Cq(hr) ) where Cq P (iir) ~ 1/Q and Cq^r) ~ Q°. 
This is simply a reorganization of the perturbative series. For instance, consider the following 
expansion of the amplitude in the theory without pions: 
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(19) 
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where A < 1/a. The series with A = and with A^O will both reproduce effective range 
theory, but differ in the location of the pole that appears at each order in the perturbative 



expansion. In the 3 Si channel, the pole of the physical amplitude is at — ip = y/MEd = 
45.7 MeV, where Ed is the binding energy of the deuteronQ. For comparison, 1/a = 36.3 MeV 
in this channel. For A = 0, the pole that appears at each order in the perturbative expansion 
will be off by 30%. For some calculations, such as processes involving the deuteron ||, a 
better behaved perturbation series is obtained by choosing 1/a + A = y/MEd- If we want to 
reproduce the expansion in Eq. (|l9|) in the theory without pions then part of Cq(ixr) must 
be treated perturbatively. 

In the theory with pions and a non- vanishing Cq(hr), the amplitude can be obtained 
from Eqs. (|i0H T3|) by substituting Co(//r) — > Cq P ({1r) and 

- A™ - [A™)'Wjjd_ ■ (20) 

In the OS scheme, the renormalization group equation makes Cq(hr) /[Cq P (ijlr)\ 2 equal to 
a constant. Therefore, Cq(hr) can be simply absorbed into the definition of D 2 ([ir). Be- 
cause of this, the value of D 2 {[Ir) extracted from fits to NN scattering data may differ 
from the value of the renormalized coupling in the Lagrangian. In the PDS scheme, the 
renormalization group equation for C p (hr) is 



2/2 



(21) 



with solution 



+ 7775 ^,np/, . \2 = constant • I 22 ) 



Therefore, breaking Co(//r) into perturbative and nonperturbative pieces results in a mani- 
festly p,R independent amplitude in PDS|. The constraint in Eq. (15) is now fiR independent. 



Integrating out the pion gives low-energy theorems for the coefficients Vi in the effective 
range expansion 



1C 



pcot5(p) = — - + y p 2 + v 2 p 4 + v 3 p 6 + v 4 p s + . . . . (23) 



x In fact, in the 3 S\ channel the fit value of Cofm,) from Eq. flTSP gives 7 = 47.3 MeV. 
2 We would like to thank Mark Wise for pointing this out to us. 
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FIG. 3. The effective field theory and Nijmegen Partial Wave analysis || values of p cot 8 are 
compared. The solid lines use pcot(5^ + S^ 1 '), the dashed lines use Eq. (^|) with the Vi from 



Ref. [ 10 1 , and the dotted lines use the values of V{ from the low-energy theorems. 



The Vi can be predicted in terms of one parameter, Cq P (/zh), which is fixed in Ref. ]T0| by 
the condition 4tt /[MCq P {^r)] + /Ir— 1/a. Corrections to these predictions are expected to 
be 30 — 50% due to higher order Q/A terms. The Vi extracted from the phase shift data 
disagree with the low-energy theorems by factors of order 5. In Fig. |3|, we see that 



the agreement of pcot(S^ + 5^) (solid lines)^ with the Nijmegen partial wave analysis is 
comparable to that of the effective range expansion with the from the fits in Refs. ||10|]TT|| 
(dashed lines). Note that our fit is more accurate at low momentum than the global fit in 
Ref. ||. However, keeping only the first five terms from the low-energy theorems (dotted 
lines) gives larger disagreement at 70 MeV. This is not surprising since the pion introduces a 
cut at p — so the radius of convergence of the series expansion of pcot (5) in Eq ( |23| ) 

is ~ 70 MeV. At p = 70 MeV, one expects large corrections from the next term in the 
series. However, the fit values of Uj give good agreement with the data even at 70 MeV. It 
is possible that uncertainty from higher order terms in the Taylor series has been absorbed 
into t>2, t>3, and t>4 in the process of performing the fits. For this reason, the uncertainty in 
the values of Vi that were found from fitting to the data may be considerable. 



3 Note that when expanded in Q, pcot5 = ip + Att /[MA { -^] -A7rA^/[M(A^) 2 ] + C(Q 3 ), 
which differs from pcot(S^ + 5^) by terms of order Q 3 . The latter expression is used since 
the parameters in Eq. (18) were fit using Eq. dl7|) . 
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To get an idea of the error in V2, we will specialize to the 3 Si channel. The Nijmegen 
phase shift analysis |12| lists two data points for p < 70 MeV: p = 21.67 MeV where 6^ SlS> = 
147.747 ± 0.010 °, and p = 48.45 MeV, where = 118.178 ± 0.021 °. Using a = 5.420 ± 



0.001 fm and ro = 1.753 ± 0.002 fm |TT[ in the effective range expansion and fitting to 
the lowest momentum data point, we find t>2 = —0.50 ± 0.52 fm 3 , where the error in a, 
ro, and pcotS have been added in quadrature. This differs by one sigma from both the 
value predicted by the low-energy theorem, v^ m = — 0.95fm 3 , and the value from the fit, 
ff 1 = 0.04 fm 3 . Since the range of the pure nucleon theory is 70 MeV, there will also be a 
~ 0.1 fm 3 error in this extraction from v% and higher coefficients. This error was estimated 
by comparing the theoretical expression for p cot 5 with the first three terms in its series 
expansion. If we instead use the higher momentum point we find t>2 = 0.03 ± 0.04 fm 3 with 
~ 0.5 fm 3 theoretical uncertainty. The uncertainty in these values of Vi is too large to make 
a definitive test of the low-energy theorems. 

In this paper, we have addressed the issue of fiR sensitivity in perturbative treatments of 
the pion. Amplitudes are \ir independent in the OS scheme, and in the PDS scheme, if part 
of Co(/j,ji) is treated perturbatively. Fits to NN scattering data were done which agree well 
at low momentum. Errors at high momentum are consistent with uncertainty from higher 
order terms if the range of the effective field theory is > 300 MeV. We conclude that there 
is no obstruction to using perturbative pions for momenta p > m n . In a future publication 
0, we will describe the renormalization procedure and OS scheme in greater detail. The 
range of the theory with perturbative pions will also be investigated. 

We would like to thank Mark Wise for many useful conversations. We also would like to 
thank H. Davioudiasl, S. Fleming, U. van Kolck, Z. Ligeti, S. Ouellette and K. Scaldefferri, 
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was supported in part by the Department of Energy under grant number DE-FG03-92-ER 
40701. 



11 



REFERENCES 



[1] S. Weinberg, Phys. Lett. B251 (1990) 288; Nucl. Phys. B363 (1991) 3; C. Ordonez 
and U. van Kolck, Phys. Lett. B291 (1992) 459; C. Ordonez, L. Ray and U. van Kolck, 
Phys. Rev. Lett. 72 (1994) 1982; Phys. Rev. C53 (1996) 2086; U. van Kolck, Phys. Rev. 
C49 (1994) 2932; D.B. Kaplan, M.J. Savage, and M.B. Wise, Nucl. Phys. B478 (1996) 
629; G.P Lepage, frucl-th/ 97060291 J.V. Steele and R.J. Furnstahl, |nucl-th/9802069l 

[2] D.B. Kaplan, M.J. Savage, and M.B. Wise, Phys. Lett. B424 (1998) 396; D.B. Kaplan, 
M.J. Savage, and M.B. Wise, |nucl-th/9802075| . 



[3] 
[4] 
[5] 
[6] 
[7] 
[8 
[9 
[10 

[11 



U. van Kolck, |nucl-th/9808007 . 

J. Gegelia, |nucl-th/9802"038; . 

J. Gegelia, |nucl-th/9806"02S . 

J.V. Steele and R.J. Furnstahl, |nucl-th/ 9808022 . 

T. Mehen and I. Stewart, CALT-68-2194. 

V.G.J. Stoks, et.al., Phys. Rev. C49 (1994) 2950. (cf . |http: / / nn-online.sci.kun.nl/N N /| ) 
D.B. Kaplan, M.J. Savage, and M.B. Wise, |nucl-th/980403l . 
T.D. Cohen and J.M. Hansen, |nucl-th/9808U3"8] . 
V.G.J. Stoks, et.al. 
th/9509032 . 



Int. Symp. on the Deuteron, Dubna, Russia, Jul 4-7, 1995. [nucl- 



[12] V.G.J. Stoks, et.al, Phys. Rev. C48 (1993) 792. 



12 



